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Abstract 

The two-terminal conductance of a random flux model defined on a square lattice is investigated numerically at the band 
center using a transfer matrix method. Due to the chiral symmetry, there exists a critical point where the ensemble averaged 
mean conductance is scale independent. We also study the conductance distribution function which depends on the boundary 
conditions and on the number of lattice sites being even or odd. We derive a critical exponent u = 0.42 ± 0.05 for square 
samples of even width using one-parameter scaling of the conductance. This result could not be obtained previously from 
the divergence of the localization length in quasi-one-dimcnsional systems due to pronounced finite-size effects. 
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1. Introduction 

The electronic transport properties of disordered 
materials depend essentially on the physical symme- 
try of the system under examination. The well known 
standard symmetry classes — orthogonal, unitary, and 
symplectic — are appropriate for describing situations 
that occur in the presence or absence of time rever- 
sal symmetry, and broken spin rotational symmetry, 
respectively. Recently, an additional chiral symmetry 
has been attracted considerable interest for disor- 
dered two-dimensional systems [1,2,3,4]. This special 
symmetry may appear if the underlying lattice is bi- 
partite [5], at which eigenfunctions can occupy only 
one of the two sub-lattices. In two-dimensional lattice 
models with diagonal disorder, the chiral symmetry is 
not observed so that for orthogonal and unitary sym- 
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metry all electronic states are localized. However, for 
non-diagonal disorder situations, like the motion of 
quantum particles subject to random magnetic fluxes 
as investigated in the present contribution, the eigen- 
values occur in pairs Te^, a characteristic feature of 
the chiral symmetry [5]. Appertaining to the chiral 
symmetry is a critical point in the center of the tight- 
binding band with a diverging localization length at 
energy E = [6,7] and a scale independent critical 
conductance. 



2. Model and Method 

We consider a two-dimensional tight-binding lat- 
tice model where the perpendicular random mag- 
netic field is introduced via complex hopping terms. 
These are chosen such that the magnetic flux per 
plaquette is given by the sum of the random Peierls 
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phases along the two bonds in the z-direction 2n(j) m = 
a m ,m+az — OL m +a x , m +a z - The corresponding Hamilto- 
nian with nearest neighbor hopping is 



H ^ ] { C rn + a x C m "T" c m ~a x C rn + 



(1) 



where the width L (^-direction) and the length L z (z- 
dircction) of the sample is given in units of the lattice 
constant a. The and c m create or annihilate a Fermi 
particle at the site m, respectively. The random fluxes 
are distributed uniformly, — f/2 < 4> m < f/"2, with 
probability density p{<f> m ) = 1//, and the randomness 
is maximal for / = 1 used in the following. The average 
magnetic flux through the system is zero. 

We calculate the dimensionless two-terminal con- 
ductance [8] for systems of length L z via 



g = Tr{iM = 



costi 2 ^/^) 



(2) 



Here, t is the transmission matrix with the Xi parame- 
terizing its eigenvalues. N is the number of open chan- 
nels. Dirichlet boundary conditions are imposed in the 
transversal direction. 



3. Conductance 

Fig. 1 confirms the theoretical prediction [3] that for 
quasi one-dimensional systems (Q1D), (Ing) decreases 
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Fig. 1. The dependence of Ing on the square root of the 
aspect ratio L z /L (length to width). 
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Fig. 2. The probability distribution p(\xi |) for odd (L = 65) 
and even (L = 66) sample widths. 

much slower for L being odd than for the same system 
with L even. This is easy to understand from Eq. (2) 
and from the mean values of the parameters Xi. In sys- 
tems with chiral symmetry, the xt assume both pos- 
itive and negative values. When ordered in absolute 
values, xi | < \x2 < . . . , it follows from symmetry rea- 
sons that (xi) = for L odd, but (a;i) = — (#2) 7^ 
when L is even. This is confirmed in Fig. 2 which shows 
the probability distribution p(|a;i|) both for L odd and 
even. For L even, p(|a;i|) converges to the Gaussian dis- 
tribution with non-zero mean when the length L z of 
thesystem increases, butp(|a;i|) remains half-Gaussian 
for L odd. 

The different form of the distribution of x\ influences 
the transmission through very long Q1D systems and 
thus the conductance. Fig. 3 shows the distribution of 
the logarithm of the conductance p(lng). In very long 
Q1D systems, the conductance is determined by the 
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O 65x1600: <ln g>=-6.205 var =29.7 
O 66x1600: <lng>=-24.4 var=49.2 




Fig. 3. The distribution of Ing. The solid lines are given by 
the function (4). Note the divergence of p(ln g) for Ing — * 
when L is odd. 
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first channel, g = cosh 2 (xi), so that 



dxip(xi) 5(\n g — 2 In cosh x\) 



(3) 



: coth(i/2)exp(-i 2 /2cr), (4) 



p(\ng) -. 
We find 



where £ = 2\n[e~ y/2 + Ve~ y - 1], and a = (if). For 
L even, p(xi) is Gaussian with a non-zero mean value. 
Consequently, p(ln g) is Gaussian, too. The singularity 
of coth(s) in Eq. (4) is eliminated by the exponential 
decrease of p{x\). However, for L odd, p(xi) possesses 
a maximum at x% = so that p(ln g) is singular for 
In g — > 0. Since this singularity survives for any length 
of the system, there is a non-zero probability to obtain 
a sample with conductance g of order of unity. This ex- 
plains the absence of exponential localization in these 
systems, we find (g) ~ L z for L odd instead as 
shown in Fig. 1. 



4. Scaling of the conductance 

Single parameter scaling theory requires that the 
conductance is a function of only one parameter in the 
vicinity of the critical point, 



(g{E,L))=F\L/i{E)]. 



(5) 



Here f(-E) is a correlation length which diverges as 
£{E) oc In Ref. [9] we have analyzed the scal- 

ing of the smallest Lyapunov exponent of Q1D systems 
with odd width and found the critical exponent v w 
0.35 ± 0.03. As also shown in Ref. [9], no scaling analy- 
sis of the Lyapunov exponent is possible for L even be- 
cause of strange finite size effects. Therefore, we verify 
in this paper the one parameter scaling for 2D square 
samples L x L with even L using the mean conduc- 
tance (Fig. 4). First, we subtract from the conductance 
the finite size correction, given by the irrelevant scaling 
term c/L, 



{g(E = 0,L))=g c -e/L, 



(6) 



where c = 2.06 and g c = 1.489 ± 0.001. Then, the data 
scale to the universal curve, (g) — c/L — F ^EL 1//y j 
with a critical exponent v — 0.42±0.05, consistent with 
our recent result for the divergence of the localization 
length in samples having odd width [9] . 



1.5 



1.45 



" 1.4 



A 
tit 
V 



1.35 

1.3 
0. 




o 


L = 


64 


■ 


L = 


92 


o 


L = 


128 


▲ 


L = 


192 


V 


L = 


256 


• 


L = 


384 


> 


L = 


512 


< 


L = 


560 



0.01 
1/L 

: 1.489 + 2.06/L 



01 



100 



10000 



EL 



i/v 



Fig. 4. The scaling of the conductance in the vicinity of 
the band center. The inset shows the finite size correc- 
tions at the critical point. After subtracting the irrelevant 
scaling term 2.06/L from the conductance, all data with 
|-E| < 0.0005 scale to a universal curve. 

5. Conclusion 



The two-terminal conductance of a random flux 
model has been investigated numerically. It shows 
critical behavior at the band center due to the chiral 
unitary symmetry of the system. The critical expo- 
nent obtained from a finite size scaling analysis agrees 
with the one that governs the divergence of the lo- 
calization length. These results support the notion of 
universality at the chiral critical point. 
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